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ABSTRACT 



1.^ ' We laser-cool single beryllium ions in a Paul trap to the ground {n = 0) quantum harmonic oscillator state with 
QJ greater than 90% probability. From this starting point, we can put the atom into various quantum states of motion 
'"^ ■ by application of optical and rf electric fields. Some of these states resemble classical states (the coherent states), 
0^ ,; while others are intrinsically quantum, such as number states or squeezed states. We have created entangled position 
and spin superposition states (Schrodinger cat states), where the atom's spatial wavefunction is split into two widely 
separated wave packets. We have developed methods to reconstruct the density matrices and Wigner functions 
^ ' of arbitrary motional quantum states. These methods should make it possible to study decoherence of quantum 
00 . superposition states and the transition from quantum to classical behavior. Calculations of the decoherence of 
CO ■ superpositions of coherent states are presented. 

O^ ' Keywords: quantum state generation, quantum state tomography, laser cooling, ion storage, quantum computation 
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1. INTRODUCTION 

t^. In a series of studies we have prepared single, trapped, ^Be+ ions in various quantum harmonic oscillator states and 
■^-^ performed measuremen±Ki)ii^hose states. In this article, we summarize some of the results. Further details are given 
S . in the original reports. EI'dBqIj 

^ . The quantum states of the simple harmonic oscillator have been studied since the earliest days of quantum 

^ ' mechanics. For example, the harmonic oscillator was among the first applications of the matrix mechanics of Heisen- 
^ , berg3 and the wave mechanics of Schrodinger. Q The theoretical interest in harmonic oscillators is partly due to the 
fact that harmonic oscillator problems often have exact solutions. In addition, physical systems, such as vibrating 
molecules, mechanical resonators, or modes of the electromagnetic field, can be modeled as harmonic oscillators, so 
C^ ■ that the theoretical results can be compared to experiments. 

A single ion in a Paul trap can be described effectively as a simple harmonic oscillator, even though the Hamil- 
tonian is actually time-dependent, so no stationary states exist. For practical purposes, the system can be treated 
as if the Hamiltonian were that of an ordinary, time-independent harmonic oscillator (see, e.g., Refs. §|,^. 

2. SYSTEM AND EFFECTIVE HAMILTONIAN 

The effective Hamiltonian for the center-of-mass secular motion is that of an anisotropic three-dimensional harmonic 
oscillator. If we choose an interaction Hamiltonian which affects only the x-motion, then we can deal with the 
one-dimensional harmonic oscillator Hamiltonian, 

Hx = hujxalax, (1) 
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where lOx is the secular frequency for the x-motion, and aj. and a^ are the creation and annihilation operators for the 
quanta of the x-oscillation mode. The constant term hujx/l, which results from the usual quantization procedure, 
has been left out for convenience. The eigenstates oi H^ are \nx), where 



H^ 



nT-hujrr 



(2) 



The internal states of the ^Be'*' ion which are the most important for the experiments are shown in Fig. n^. We are 
mostly concerned with the hyperfine-Zeeman sublevels of the ground 2s '^Si/2 electronic state, which are denoted by 
1 2s ^ 5*1/2 ; F, Mp) , where F is the total angular momentum, and Mp is the eigenvalue of F^ . The ^Be nucleus has spin 
3/2. Of chief importance are |2s ^51/2; 2, 2), abbreviated as 1 1), and |2s ^6*1/2; 1, 1), abbreviated as 1 1)- The hyperfine- 
Zeeman sublevels of the 2p ^Pj (J = 1/2 or 3/2) fine-structure multiplet also play a role, either as intermediate states 
in stimulated Raman transitions or as the final states in resonantly-driven single-photon transitions used for laser 
cooling or state detection. We denote these states by \2p^Pj; F, Mp)- The energy separation of 1 1) and 1 1) is Tiloq, 
where wo ~ 27r x 1.250 GHz. They are coupled by laser beams Rl and R2, through the intermediate \2p^Pi/2\ 2, 2) 
state. The frequency detuning from the intermediate state is A, where A w — 27r x 12 GHz. 

If we restrict the internal states to the space spanned by | ]) and | [), then the internal Hamiltonian can be written 

hujQ 



as 



H\i - 



(3) 



where a^ is a Pauli spin matrix whose nonzero matrix elements are(t \(Tz\ T) = +1 ^^^d (| jcr^l [) = —1. 

The total effective Hamiltonian for the system consisting of the 
internal | f) and | J,) states and the x motional degree of freedom is 

H = Hn + H, + Hi^t, (4) 

where Hmt is the effective interaction Hamiltonian coupling due to 
the two laser beams Rl and R2 in Fig. ||. We define an interaction- 
picture operator O^ in terms of a Schrodinger-picture operator O^ 
as 

0\t) = e'"'^^^/^'0^{t)e-''"'^i^'^. (5) 

The effective interaction Hamiltonian in the interaction picture and 
the rotating wave approximation is 



'1/2 




-ff;„t = 



hg [a+e 



iv('ll:+0-o:)—iSt 



^ — iri(al.+a^)+iSt 



(6) 



Figure 1. Internal and motional energy 
levels of a ^Be"*" ion. The | t) and | [) states 
are particular hyperfine-Zeeman components 
of the ground ^5*1/2 state, separated in energy 
by Tildq. Each internal state can exist in any 
of a ladder of vibrational energy states jrij.), 
where n^, = 0, 1, 2, . . ., separated in energy by 
huj^, where lu^ k,2tt y. 11.2 MHz. 



where g is the interaction strength, 5 is the detuning of the fre- 
quency difference of the two laser beams with respect to o^o, and 
77 = k^Jh/ {2mu!x) is the Lamb-Dicke parameter, where k is the 
magnitude of the difference between the wavevectors of the two laser 
beams, and m is the mass of the ion. The nonzero matrix elements 
of (7+ and a^ are (t |cr+ 1 J,) = (J, |o-_ 1 1) = 1. 

The detuning 5 can be tuned to multiples oi oJxj 5 = {n' — n)ujx, 
so as to resonantly drive transitions between 1 1, n) and 1 1, n'). We 
refer to the 6 = resonance as the carrier, the d = —lUx resonance as 
the first red sideband, and the 5 = +ujx resonance as the first blue 
sideband. 



The signal that is detected in the experiments is the probability 
P|(i) that the ion is in the | |) internal state after a particular 
preparation. If we irradiate the ion at time t with circularly polarized light, resonant with the electronic transition 
from I i) to |2p^F3/2; 3, 3), there will be a high fluorescence intensity if the 1 1) state is occupied, since the selection 
rules only allow that upper state to decay back to the 1 1) state, so it can continue to scatter photons. This transition 
is called a cycling transition. If the ion is in the 1 1) state when it is irradiated with the same light, it will scatter a 
negligible number of photons. Thus, if we repeatedly prepare the ion in the same way, apply radiation resonant with 
the cycling transition, and detect the fluorescence photons, the average signal will be proportional to Pi{t). 



3. CREATION AND PARTIAL MEASUREMENTS OF QUANTUM STATES 

3.1. Fock states 

The ion is prepared in the n — state by Raman cooHng.EJ (From now on we drop the x label on n.) Raman cooling 
consists of a sequence of laser pulses on the red sideband, driving | i, n) to | t, 't- — 1) transitions, followed by laser 
pulses which recycle the ion back to the ||, n— 1) state. The probability of heating due to recoil during the recycling 
step is small. At the end of the sequence, the ion is in the n — state more than 90% of the time. 

A Fock state is another name for an n-state. Higher-n Fock states are prepared from the n — state by a 
sequence of 7r-pulses on the blue sideband, red sideband, or carrier. For example, the 1 1, 2) state is prepared by using 
blue sideband, red sideband, and carrier 7r-pulses in succession, so that the ion steps through the states | i, 0), 1 1, 1), 

|i,2),and|T,2). 

If the atom is initially in the 1 1, n) state, and the first blue sideband is driven, it will oscillate between that state 
and the 11,?^+ 1) state. The probability of finding it in the | i, n) state at a time t is 



nW = I [1 + cos(2rj„,„+it)e-''"*] , 



(7) 



where fi 



n,n+l 1 



the Rabi flopping rate, is a function of the laser intensities and detunings and n, and 7„ is a damping 
factor which is determined empirically. Thus, the frequency of the oscillations of P^ (t) is a signature of the initial 
value of n. Figure g(a) shows Pi{t) for an initial | i, 0) state. Figure ||(b) shows the observed ratios of Rabi frequencies 
compared with the values calculated from the matrix elements of the interaction Hamiltonian [Eq. m)]. 
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Figure 2. (a) P|(i) for an initial | i,0) number state driven by the first blue sideband. The solid line is a fit to 
an exponentially damped sinusoid, (b) Observed ratios of the Rabi frequencies Qn.n+i/^o,i for different values of 
the initial n. The lines represent calculated values for different values of the Lamb-Dicke parameter rj. The value of 
?7, inferred from other measured quantities, was 0.202 ± 0.005. (c) Pi{t) for a thermal state. The solid line is a fit 
of the data (dots) to a sum of Fock states having a thermal distribution. The fitted value for the mean value of n 
is n = 1.3 ± 0.1. The inset shows the amplitudes of the Fock state components (bars) with a fit to an exponential, 
corresponding to n = 1.5 ± 0.1 (line). 



3.2. Thermal states 

A thermal state of the motion of the ion is not a pure state, but rather must be described by a density matrix, 
even though there is only one ion. The statistical ensemble which is described by the density matrix is generated by 
repeatedly preparing the state and making a measurement. 



The ion is prepared in a thermal state by Doppler laser cooling.E2l The temperature of the distribution can be 
controlled by changing the detuning of the cooling laser. When the ion's state is not a Fock state, P^ (t) has the form 



nw 



1 + ^ P„ cos(2fi„,„+it)(; 



-l,^t 



(8) 



where P„ is the probability of finding the ion in the state |7i). Figure 0(c) shows Pi{t) for a thermal state. 



3.3. Coherent states 

Coherent states of the quantum harmonic oscillator were introduced by Schr6dinger,tij with the aim of describing a 
classical particle with a wavefunction. A coherent state \a) is equal to the following superposition of number states: 



a 



-hM'S-^ 






(9) 



In the Schrodinger picture, the absolute square of the wavefunction retains its shape, and its center follows the 
trajectory of a classical particle in a harmonic well. The mean value of n is rl = |ap. A coherent state can be created 
from the n = state (a special case of coherent state) by applying a spatially uniform classical force (see Appendix 
A). The drive is most effective when its frequency is uj^. Another method is to apply a "moving standing wave," that 
is, two laser beamis_differing in frequency by uj^ and differing in propagation direction, so that an oscillating dipole 
force is generated.E3 We have used both of these methods to prepare coherent states. Figure y(a) shows Pi{t) for a 
coherent state. This trace exhibits the phenomenon of collapse and revival.E3 
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Figure 3. (a) Pi{t) for a coherent state. The solid line is a fit of the data (dots) to a sum of Fock states having 
a coherent state distribution. The fitted value for n is 3.1 ± 0.1. The inset shows the amplitudes of the Fock state 
components (bars) with a fit to a Poissonian distribution, corresponding to n = 2.9 ± 0.1 (line), (b) Pi{t) for a 
squeezed state. The solid line is a fit of the data (dots) to a sum of Fock states having a squeezed-state distribution. 
The fitted value for /3 is 40 ± 10, which corresponds to n « 7.1 



3.4. Squeezed states 

A "vacuum squeezed state'Lcan be created by applying an electric field gradient having a frequency of 2uJx to an ion 
initially in the n = state.tj Here, the ion was irradiated with two laser beams which differed in frequency by 2a;a;. 
This has the same effect. The squeeze parameter /3 is defined as the factor by which the variance of the squeezed 
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Figure 4. Creation of a Schrodinger cat state. See text for details. 

quadrature is decreased. It increases exponentially with the time the driving force is applied. The probability 
distribution of n-states is nonzero only for even n, for which, 

(2n)!(tanhr)^" 
^'" - ^ (2^^ ' ^^^> 

where /3 = e^^, and iV is a normalization constant. Figure ra(b) shows Pi{t) for a squeezed state having /? = 40± 10. 

3.5. Schrodinger cat states 

The term "Schrodinger cat" is taken here to denote an entangled state which consists of two coherent states of motion 
correlated with different internal atomic states. A simple example is 

I*) = mi!>±lk^. (u) 

Figure |j shows how such a state is created, (a) The ion is prepared in the | J,, 0) state by Raman cooling, (b) A 
7r/2-pulse on the carrier generates an equal superposition of 1 1, 0) and | Ti 0). (c) The displacement beams generate 
a force Fi, which excites the component in the | f) internal state to a coherent state \a). Due to the polarizations 
of the displacement beams, they do not affect an atom in the | J,) state, (d) A 7r-pulse on the carrier exchanges 
the 1 1) and 1 1) components, (e) The displacement beams generate a force F2, which excites the component in the 
I t) internal state to a coherent state jae*"^), where the phase (j) is controlled by an rf oscillator. The state here is 
analogous to Schrodinger's cat. (f) The 1 1) and | "f) components are combined by a 7r/2-pulse on the carrier. At this 
point, the radiation is applied on the cycling transition, and the signal is recorded. 

The predicted signal, for a particular value of (j), is 

P^(0) = ^ [1 - ce-"'(i-™''"*) cos(a2 sin0)] , (12) 

where c = 1 in the absence of decoherence. Figure shows experimental data and fits to Eq. (|l2) for various values 
of a. In (d), c is clearly less than one, indicating decoherence, although the source is not yet determined. 

4. COMPLETE MEASUREMENTS OF QUANTUM STATES 

The measurements described in the previous Section determine only the n-state populations or probabilities, and 
therefore do not provide a complete description of the motional states. The density matrix p does provide a complete 
description of a state, whether it is a pure or a mixed state. The Wigner function PF(a) also provides a complete 
description. The Wigner function resembles a classical joint probability distribution for position and momentum in 
some cases. However, it can be negative, unlike a true probability distribution. We have demonstrated experimental 
methods for ^constructing the density matrix or the Wigner function of a quantum state of motion of a harmonically 
bound atom.cl 
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Figure 5. Measured and fitted Pi{<t>) interference signals for tlie Schrodinger cat state. The curves in (a) through 
(d) represent measurements for different values of the time during which the displacement beams excite the coherent 
states. The solid lines are fits to Eq. (I2h, allowing a to vary. In (d), a term representing loss of contrast was included 
in the fit. The curve in (d) represents a superposition of two xq ~ 7 nm wave packets with a maximum separation 
of 4axo ~ 80 nm. Curve (e) represents a calculation for a pair of coherent states with a — Q. 



Both of these methods depend on controUably displacing the state in phase space, applying radiation to drive 
the first blue sideband for time t, and then measuring P^. The averaged, normalized, signal is 



Piit,a) = ^ 



1 + ^ Qk{a) cos{2Vtk,k+it)t 



-ikt 



k=0 



(13) 



where the complex number a represents the amplitude and phase of the displacement, and Q{a) is the occupation 
probability of the vibrational state |fc) for the displaced state. 

If the Qkia) coefficients are derived for a series of values of a lying in a circle, 



ifp = |a|exp[i(7r/7V)p], 



(14) 



where p 



-N, 



,N — 1^ then the density matrix elements pnm can be determined for values of n and m up to 



iV — 1. The details of the numerical procedure are given in Ref. ^. 

Figure @(a) shows the reconstructed density matrix amplitudes for an approximate n = 1 state. Figure ra(b) 
shows the reconstructed density matrix for a coherent state having an amplitude \f3\ w 0.67. 

The Wigner function for a given value of the complex parameter a can be determined from the sun£3ll3'0'E3 
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Figure 6. (a) Reconstructed number-state density matrix amplitudes \pnm\ for an approximate n = 1 nmnber 
state, (b) Reconstructed amplitudes \pnrn\ and phases Q{pnrn) of a coherent state. 

Figure H shows the reconstructed Wigner function for an approximate n = 1 state. The fact that it is negative in 
a region around the origin highlights the fact that is a nonclassical state. Figure H shows the reconstructed Wigner 
function for a coherent state with amplitude |/3| w 1.5. It is positive, which is not surprising, since the coherent state 
is the quantum state that most closely approximates a classical state. 
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Figure 7. Surface and contour plots of the reconstructed Wigner function M^(a) of an approximate n = 1 number 
state. The negative values of W{a) around the origin highlight the nonclassical nature of this state. 




Figure 8. Surface and contour plots of the reconstructed Wigner function 14/^(0;) of a coherent state. The width is 
in good agreement with the expected minimum-uncertainty value. 



5. DECOHERENCE OF A SUPERPOSITION OF COHERENT STATES 

Decoherence, the decay of superposition states into mixtures, is of interest because it leads to the familiar classi- 
cal physics of macroscopic objects (see, e.g. Ref. |l9|) . Experiments are beginning to probe the regime in which 
decoherence is observable, but not so fast that the dynamics are classical. One example is the observation of the 
decoherence of a superposition of mesoscopic (few-photon) quantum states of an electromagnetic field having distinct 
phases by Brune et alS3 In this experiment, as in the Schrodinger cat experiment of Ref. y, the mesoscopic system 
was entangled with an internal two-level system of an atom. 

A simpler system, which has been well-studied theoretically by others, EilESE^O is a superposition of two coherent 
states of a harmonic oscillator, without entanglement with a two-level system. The system interacts in some way 
with the environment, leading to decoherence of the superposition state. 

Here we carry out a simple calculation of the decoherence of a superposition of two coherent states of a harmonic 
oscillator. At time t — 0, the system is in an equal superposition of two coherent states. 



l*(0)) 



1 (^^em 
V2 V 



«i(0)) 



p»e2(o) 



«2(0)) 



(16) 



where the notation is the same as in Appendix A, and the interaction picture is used. (It is assumed that the two 
components of the state vector are nearly orthogonal. Otherwise the normalization factor is different.) The system 
is subjected to a random force, which is uniform over the spatial extent of the system, for example, a uniform electric 
field acting on a charged particle. In Appendix A, wc show that a single coherent state, subjected to such a force, 
will remain in a coherent state. It follows from the linearity of the Schrodinger equation that the state described by 
Eq. (|l^) will, at a later time, still be in an equal superposition of coherent states, although ai, a2, 9i, and 02 will 
change. Each time the system of Eq. (Hq) is prepared, ai, a2, 9i, and 02 will have a different time dependence, due 
to the random force. We can study the decoherence as a function of time by carrying out an ensemble average. 

Although both ai and a2 change with time, their difference does not, since they each change by the same amount. 
This can be seen from Eq. (Eq). We define Aa = ai — a2 = constant. 



Decoherence is due to the changes in the phase difference [6*2 (t) — ^i(i)]- We define the change in the phase 
difference to be 

Aeit) = [02 (t) - 01 (t)] - [92 (0) - e, (0)] . (17) 

Then, from Eq. iM), 



mt) = I 






(18) 



Aa / /(t')e"'"^* dt' + Aa* / /(i')e"^==* dt 
Jo Jo 

The double-integral term in Eq. (El) does not contribute to A9(t), since it does not depend on either a(0) or 9{0) 
so it makes the same contribution to 6*2 (i) and to 9i{t). The square of A9{t) is 



imt)f = \ 



(Aa)^ / / f{t')f{t")e''^-' -"^-* dt'dt" + 2|Aa|2 



"'0 
t rt 



"'0 



^'•'f'^"^^i'^-t'-''^-t" dt'dt" 



f{t')f{t")e^' 



+ {Aa*f f f f{t')f{t") 
Jo Jo 



'\ f(i"\gi^^t'+iu:^t"^^/^^/l 



dt'dt" 



(19) 



and its ensemble average is 



iA9it)y 



t l<t 

2 / / / -cu'\ -cW'W „-i^xt' -iuj.j:t''u'i4-" 



{AaY I / {f{t')f{t"))e 
Jo 

■t rt 
2 



t i<t 

2 / / I f I -^l\ f { -l-l'W ^i^:j:t' —iijJ^t" 111 Jill 



dt'dt" + 2\Aa\-' I / {f{t')f{t"))e 
Jo 



dt'dt" 



(Aa* 



^0 



if it')f{t"))e"^^*'+''^-'" dt'dt" 



(20) 



where the angular brackets denote an ensemble average. Thus, ( {A9{t)) ) depends on the autocorrelation function 

of the random function /(i), which is proportional to the force. If we assume that / is a stationary random variable, 
then its autocorrelation function 



R{T) = {f{t)f{t + r)) 



(21) 



exists and is independent of t. According to the Wiener-Khinchinc theorem, the power spectrum w^v), defined for 
1/ > 0, is 

/•oo 

w{v) == 4 / R{t) cos{2T:vT)dT. (22) 

A white-noise spectrum corresponds to R{t) cx 5{t), where 6{t) is the Dirac delta function. If we let R{t) = CS{t), 
then 



{A0{t)y 



c 



2 / / xuii ^'^^-«'^^*'-»'^=^*"di'di" + 2|AaP 



(AaY I / Sit"-t')e 
Jo 
■t rt 



^0 



^11 ^i\„iuj^t'~iu^t"^^ij^^i, 



5{t" - t')e" 



{Aa*f f f 5{t" -t')e"^'^''+''^'''" dt'dt" 
Jo Jo 



C 



(Aa)- / e 



f e-2»'^==*'dt' + 2|Aa|2 f dt' + {Aa*f f e 
^0 Jo Jo 



^ \2 I ^2iuJxt jj-l 



(23) 



The first and third integrals on the right-hand-side of Eq. (^3|) oscillate, but are bounded, while the second one grows 
with time, so, for ujxt ^ 1, 

{{A9{t)f) « lc\Aa\H. (24) 

/ 2\^/^ 

Hence, the decoherence time, that is, the time required for the rms phase difference ({A9{t)) ) to increase to 

about 1 radian, is on the order of 2/ (C|Aap). 

Using the same approximations, we can calculate the rate of change of the amplitude ai{t) for a single coherent 
state, where i is 1 or 2. From Eq. (Esh in Appendix A, 



a,(t)-a.(0) = z / /(t')e'"^*'di'. 



(25) 



and 

|a,(t)-«»(0)p- / [ f{t')f{t")e'^-''-'^-'"dt'dt". (26) 

JO JQ 

If f{t) has a white-noise spectrum, as considered previously, the ensemble average is 

{\a,{t) ~ a,(0)p) =C f f Sit" - t')e'"-''-'^-'"dt'dt" = C f dt' ^ Ct. (27) 

Jo Jo JQ 

In the time, t « 2/(C|Aq;P), required for decoherence, the rms change in ai is V2/|Aq;|. Consider the case 
ai(0) — — a2(0). The energy of a single coherent state is proportional to the square of its amplitude (oc |aip). The 
fractional change in the energy is 



Ala,: 



2|a,(t)-a,(0)| ^ 2^/2 ^ V2 

a,(0) ^ |Aa||a,(0)| |«,(0)|2' ^ ' 

which becomes small for |Q;i(0)p 3> 1. Thus, we see that coherent superpositions of macroscopic (|ai(0)p ^ 1) states 
decohere much more quickly than they change in energy. 

We can make a more quantitative statement about the decoherence rate by considering the decay of the off- 
diagonal density matrix elements. For an initial pure state 2~^''^{\ai) + |a2)), the density matrix is 

P=2 (l"i)("il + laiXaal + |a2)(ai| + |a2)(a2|) . (29) 

If Aa ^ 1, so that the phase changes are more important than the changes in ai and a2, then the off-diagonal matrix 
element |a2)(ai| evolves to e'^^'^*-'|a2)(ai|. Let us assume that the random variable A9 at time t has a Gaussian 
distribution P{A9). Given the variance of A9 from Eq. (p4|), P{A9) must have the form 

We evaluate the ensemble average of e*'^^'*^ 

^Aeit)\^ r F(A0)e'^«Wd(A0) = exp f-^^^^V (31) 



This yields the basic result, previously derived by others,EifE3E3cil that the off-diagonal matrix element decays 
exponentially in time, with a time constant that is inversely proportional to |Aap. 

An extension of the experiment of Ref. H, with a random electric field deliberately applied, might be used to 
verify these calculations.B A means of simulating thermal, zero-tempemture, squeezed, and other reservoirs by 
various combinations of optical fields has been discussed by Poyatos et al£B 
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APPENDIX A. THE FORCED HARMONIC OSCILLATOR 

If a coherent state is subjected to a spatially uniform force, it remains a coherent state, though its amplitude 
changes. This exact resul 



^_,_ which is independent of the strength or time-dependence of the force, seems to have 

been discovered by HusimiEj and independently by KernerH3 For other references on the forced quantum harmonic 
oscillator, see Ref. g^ The qualitative result stated above is enough to establish that applying a force to the n = 
state will generate a finite-amplitude coherent state, which was required in Sec. ^J5. However, in order to study the 
decoherence of a superposition of two coherent states, it is useful to have an exact expression for the time-dependent 
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state vector, given that the state vector is an arbitrary coherent state at t = 0. Since the pubhshed results of which 
we are aware do not give this particular result, we give here an elementary derivation. 

We consider a Hamiltonian, H — Hq + V(t), where Hq = fiLOxol.ax is the Hamiltonian of a one-dimensional 
harmonic oscillator, and V{t) = —xFit) is a time-dependent potential. Here, F{t) is a real c-number function of 
time and corresponds to a spatially uniform force, and x ~ XQ{ax + aj,), where xq = ('h/2mujx)^''^ . For example, if 
the particle has charge q, and a uniform electric field xE^it) is applied, then F[t) = qEx{t). 

It is convenient to switch to the interaction picture, where an interaction-picture state vector |^(i))^ is related 
to the corresponding Schrodinger-picture state vector |^(t))"^ by 

|*(t))-f = e'-f^°*/''|*(t))^, (32) 

and an interaction- picture operator O^ {t) is related to the corresponding Schrodinger-picture operator 0^{t) by 

The equation of motion for an interaction-picture operator is 

lo'it)=-^[0\t),Ho], (34) 

so 

j^aiit) = -^[ai{t),Ho]^-iu;xa'it), (35) 

aiit) = ai(0)e— ^*, (36) 

al\t) = at'(0)e-^*. (37) 

For simplicity, we let a = a;J(0) = af and a^ = aj. (0) — aj, . The Schrodinger equation in the interaction picture is 

||vl>(t))/ = .^v'imit))' = »/(i) (ae--^* + ate-^*) |vl/(i))^ (38) 

where f{t) = xoF{t)/h. Note that if f{t) = 0, |*(t))^ is constant. 

We impose the initial condition that |^(0))^ is a coherent state. We make the assumption that |^(t))^ is also a 
coherent state for i > 0, 

\^{t)y = e^'^'^\a{t)) EE e^m^-h'^m^ V ^\n), (39) 

n=o ^^- 

where the \n) are the eigenstates of al.ax with eigenvalue n at t = 0, and d{t) is an arbitrary real function, so 
that ^(^(t)l^(i))^ = 1. (This assumption must be verified later by substitution of the resulting solution back into 
Eq. (p8|).) The problem reduces to finding the complex function a{t) and the real function 9{t). 

The right-hand-side of Eq. @) is 






h 



oc 



z/(t)e'^(*)e-^l"WI^ e— * ^ " ''/^ " \n - I) + e^^* ^ " '''V!^' \n + 1) (40) 



n— 1 ^ n— 

oc 



^o{t),-i\c.it)\- L-^^^t ^ ^!1M|„\ _u .-.* V^ "" '(^)^/" l 



z/(t)e^'(*)e-il"(*)l' e-^'^-* ^ ^^y^|n) + e^^^' ^ 



.^1 V(^^^l)! 
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The left-hand-side of Eq. (pi) is 
5,._,r .d9 



n— ^ 



2 fit 



1 , ^do* 



a"(t) 



2 ^ d* ; 2^^^ 



We equate the coefficients of \n) in Eqs. (|4^) and ( pl| ) and divide by a common factor to obtain 



I dt 



if{t) (e-'"-*a2(t) + e'^-'n) = i^a(f) + 



I da . , , 1 , s da* \ / X da 



which can be rearranged as 



da 



-dt-'-f^'^' 



-W? - i«W«* W$ - i"'W^ - ^a^W/(i)e--* = 



dt 2 



dt 2 



dt 



(41) 



(42) 



(43) 



In order for Eq. (Bs) to be true for aU n, the expression in the square brackets, which multipUes n, must equal zero. 
Thus, we have the first-order differential equation for a{t), 



^-^/(t)e-^* 



0, 



a{t) =a{0)+i I f{t')e"^^''dt'. 



which has the solution 



If the expression in brackets is set to zero, Eq. ( |43| ) becomes, after division of both sides by ia(t), 

d9 i ^, .da i , .da* , ^ „, ^ ,,, < 
d^+2"(^)d^+2"(*)^-"(*)-^(*)^ "^«- 



(44) 
(45) 

(46) 



The imaginary part of Eq. (46) is satisfied as long as Eq. (Bj) is, so it provides no new information. The real part 
of Eq. (^) is a first-order differential equation for 9(t): 



'di 



— f{t) [Re Q(t) cos ujxt + Im Q;(t) sin uj^t] 



Rea(O) coswjjt + ImQ;(0) sinci;a;t — cos Wajt / /(t') sin oj^^t'dt' -f sin oj^^t / f {t') cos ui^t'dt 



Rea(O) cosw^jt + Ima(O) sinwa;t + / /(t') sin [aj:r(t — t')] dt' 



(47) 



= fit) 

= fit) 
where Re and Im stand for the real and imaginary parts. Integrating Eq. (E7h, wc obtain 

0{t)^0{Q)+Rea{O) f f{t')coaLjJdt'+Ima{0) f f {t') sin lu J dt' + j dt'f{t') f dt" f{t") sin [uj^{t' ~ t")] . (48) 
Jo Jo Ja Jo 



Substitution back into the Schrodinger equation [Eq. (p8|)] verifies the solution and justifies the original assumption 
[Eq. (PI)]. 



REFERENCES 

1. C. Monroe, D. M. Meekhof, B. E. King, S. R. Jeffcrts, W. M. Itano, and D. J. Wineland, "Resolved-sideband 
Raman cooling of a bound atom to the 3D zero-point energy," Phys. Rev. Lett. 75, pp. 4011-4014, 1995. 

2. D. M. Meekhof, C. Monroe, B. E. King, W. M. Itano, and D. J. Wineland, "Generation of nonclassical motional 
states of a trapped atom," Phys. Rev. Lett. 76, pp. 1796-1799, 1996. 



12 



3. C. Monroe, D. M. Meckhof, B. E. King, and D. J. Wineland, "A 'Schrodinger cat' superposition state of an 
atom," Science 272, pp. 1131-1136, 1996. 

4. D. Leibfried, D. M. Meekhof, B. E. King, C. Monroe, W. M. Itano, and D. J. Wineland, "Experimental deter- 
mination of the motional quantum state of a trapped atom," Phys. Rev. Lett. 77, pp. 4281-4285, 1996. 

5. D. Leibfried, D. M. Meekhof, B. E. King, C. Monroe, W. M. Itano, and D. J. Wineland, "Experimental prepa- 
ration and measurements of quantum states of motion of a trapped atom," J. Modem Opt. (in press) . 

6. W. Heisenberg, "Uber quantentheoretische Umdeutung kinematischer und mechanischer Beziehungen," Z. Phys. 
33, pp. 879-893, 1925. 

7. E. Schrodinger, "Quantisierung als Eigenwertproblem II," Ann. Phys. (Leipzig) 79, pp. 489-527, 1926. 

8. P. J. Bardroff, C. Leichtle, G. Schrade, and W. P. Schleich, "Endoscopy in the Paul trap: Measurement of the 
vibratory quantum state of a single ion," Phys. Rev. Lett. 77, pp. 2198-2201, 1996. 

9. G. Schrade, V. I. Man'ko, W. P. Schleich, and R. J. Glauber, "Wigner functions in the Paul trap," Quantum 
Semiclass. Opt. 7, pp. 307-325, 1995. 

10. S. Stenholm, "The semiclassical theory of laser cooling," Rev. Mod. Phys. 58, pp. 699-739, 1986. 

11. E. Schrodinger, "Der stetige Ubergang von der Micro- zur Makromechanik," Die Naturwissenschaften 14, 
pp. 644-666, 1926. 

12. D. J. Wineland, J. C. Bergquist, J. J. Bolhnger, W. M. Itano, F. L. Moore, J. M. Gilhgan, M. G. Raizen, 
D. J. Heinzen, C. S. Weimer, and C. H. Manney, "Recent experiments on trapped ions at the National Institute 
of Standards and Technology," in Laser manipulation of atoms and ions, E. Arimondo, W. D. Phillips, and 
F. Strumia, eds., pp. 553-567, North- Holland, Amsterdam, 1992. 

13. J. H. Eberly, N. B. Narozhny, and J. J. Sanchez-Mondragon, "Periodic spontaneous collapse and revival in a 
simple quantum model," Phys. Rev. Lett. 44, pp. 1323-1326, 1980. 

14. D. J. Heinzen and D. J. Wineland, "Quantum-limited cooling and detection of radio-frequency oscillations by 
laser-cooled ions," Phys. Rev. A 42, pp. 2977-2994, 1990. 

15. A. Royer, "Measurement of the Wigner function," Phys. Rev. Lett. 55, pp. 2745-2748, 1985. 

16. H. Moya-Cessa and P. L. Knight, "Series representations of quantum-field quasiprobabilities," Phys. Rev. A 48, 
pp. 2479-2485, 1993. 

17. S. Wallentowitz and W. Vogel, "Unbalanced homodyning for quantum state measurements," Phys. Rev. A 53, 
pp. 4528-4533, 1996. 

18. K. Banaszek and K. Wodkiewicz, "Direct probing of quantum phase space by photon counting," Phys. Rev. 
Lett. 76, pp. 4344-4347, 1996. 

19. W. H. Zurek, "Decoherence and the transition from quantum to classical," Phys. Today 44 (10), pp. 35-44, 
1991. 

20. M. Brune, E. Hagley, J. Dreyer, X. Maitre, A. Maali, C. Wunderlich, J. M. Raimond, and S. Haroche, "Observing 
the progressive decoherence of the 'meter' in a quantum measurement," Phys. Rev. Lett. 77, pp. 4887-4890, 
1996. 

21. D. F. Walls and G. J. Milburn, "Effect of dissipation on quantum coherence," Phys. Rev. A 31, pp. 2403-2408, 
1985. 

22. J. P. Paz, S. Habib, and W. H. Zurek, "Reduction of the wave packet: Preferred observable and decoherence 
time scale," Phys. Rev. D 47, pp. 488-501, 1993. 

23. B. M. Garraway and P. L. Knight, "Evolution of quantum superpositions in open environments: Quantum 
trajectories, jumps, and localization in phase space," Phys. Rev. A 50, pp. 2548-2563, 1994. 

24. P. Goetsch, R. Graham, and F. Haake, "Schrodinger cat states and single runs for the damped harmonic 
oscillator," Phys. Rev. A 51, pp. 136-142, 1995. 

25. J. F. Poyatos, J. I. Cirac, and P. ZoUcr, "Quantum reservoir engineering with laser cooled trapped ions," Phys. 
Rev. Lett. 77, pp. 4728-4731, 1996. 

26. K. Husimi, "Miscellanea in elementary quantum mechanics II," Prog. Theor. Phys. 9, pp. 381-402, 1953. 

27. E. H. Kerner, "Note on the forced and damped oscillator in quantum mechanics," Can. J. Phys. 36, pp. 371-377, 
1958. 

28. Y. Nogami, "Test of the adiabatic approximation in quantum mechanics: Forced harmonic oscillator," Am. J. 
Phys. 59, pp. 64-68, 1991. 



13 



